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Abstract 

We investigate the relativistic Fermi motion effect in the case of J/tp production in 
various hadron colliders. A light-cone wave function is adopted to represent the J ftp 
final state. The change in the confinement parameter which sets a scale for the size 
of the final state, allows one to see the effect in an explicit manner. While the effect 
has considerable influence on the fragmentation probabilities and the differential 
cross sections, the total cross sections are essentially left unchanged, such a feature 
is in agreement with the momentum sum rule which the fragmentation functions 
should satisfy. 
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1 Introduction 



Evaluation of the J / ip cross section at Tevatron energies has been one of the in- 
teresting problems of QCD in theory [1,2,3] and in experiment [4]. Predictions 
using the QCD calculations where the production of J / ip is assumed to occur 
in color singlet form, fails to agree with the experimental results. To bring 
about the agreement the mechanism of color octet is introduced in which the 
bound state is produced originally in the color octet state at the production 
point [5]. Then by emitting a soft gluon the colored object is transformed into 
a color singlet state [6]. 
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Generally in the formation of quarkonia bound state it is assumed that the 
constituents are not relativistic and therefore they are let to fly together within 
the bound state ignoring their respective motion. Here we propose that the 
constituents within such states specially the charm quark is not that heavy to 
let one to ignore its relativistic effects particularly in relation to the energies 
at which such a particle is produced. Therefore the effect should be accounted 
for in the production process in a more explicit and accurate form. The results 
of such a study may shed light on the problem of J / ip production cross section. 

In this work, we introduce the Fermi motion into the J/ ip production in direct 
fragmentation process using a light-cone wave function. Since our aim is to 
show the size of the effect, we have not included other contributions. We 
demonstrate the enhancement of the fragmentation function due to this effect. 
With such a significant enhancement, we evaluate the differential cross section 
times the branching ratio and the total integrated cross section for the process 
pp — > c — > J/ip — > at the Tevatron Run I energies and compare them 

with the CDF data. We also present similar results for the case of the RHIC, 
the Tevatron Run II and the CERN LHC pp collisions. 



2 The light-cone wave function and the bound state formation via 
fragmentation 



There have been different approaches to introduce the Fermi motion into the 
production and decay processes of various meson states and quarkonia using 
different forms of wave functions [7] . The light-cone wave functions have been 
employed in the case of charmonia to study their production in B-decay [8] 
and in photoproduction [9]. Here motivated by harmonic oscillator model, we 
have picked up a wave function in the light-cone quantization to represent the 
J/ip bound state. It has the following form [7] 



ip c -*j/^(xi,x 2 ,qT) = A^j/^exp 
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where A c ^j/^ is the normalization coefficient, m is the quark (anti-quark) mass 
and qT is the transverse momentum of the constituents. The x's are the energy 
momentum ratios and finally the parameter (3 is known as the confinement 
parameter which controls the width of the wave packet representing the bound 
state. The normalization condition is 
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where 
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and 



[d 2 q T ] = J| 167r 3 5 2 [ ^ q Ti 

i=i i=i 

The sum is over all Fock states and helicities. 
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With the choice of the wave function (1), we fix all degrees of freedom (trans- 
verse and longitudinal) of the constituents within the J/ip bound state by 
matching them with those in the matrix elements relevant to the fragmenta- 
tion function and vary the parameter (3 to change the size of the wave packet 
representing the bound state. Benefiting such a method, in the leading order 
perturbative regime, the fragmentation functions for J/ip production without 
and with the Fermi motion are obtained as follows [10]: 

(a) Fermi motion off 

In this case the transverse momenta of the constituents are set equal to zero 
and that the longitudinal components are chosen to be equal. The confinement 
parameter is (3 = in this case. The fragmentation function is obtained as 
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where a s is the strong interaction coupling constant and Cp is the color factor. 
The quantity (kp 2 ) is the average transverse momentum squared of the initial 
state heavy quark, the parameter £ is defined as £ = (k T 2 )/m 2 and finally 
F(x) is given by 



F(z) 
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(b) Fermi motion on 

Here the different momentum components of the constituents, i.e. the quan- 
tities q = \qr\ and x = X\ = 1 — £2 in the matrix elements and the wave 
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function (1) are integrated over. The only remaining parameter is the con- 
finement parameter (3. The fragmentation function in this case is obtained 
as 
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The function G(z) reads 
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Here we have defined r\ = q 2 /m 2 . In our study we have set the charm quark 
mass equal to 1.25 GeV. 

The fragmentation functions (5) and (7) provide a comparison of the cases of 
the Fermi motion off and on. The two functions coincide at sufficiently low (3. 
As (3 increases, (7) raises considerably at the peak region which leads to the 
respective increase in fragmentation probability. To choose the (3 value, we 
note that different values within the range of (3 =0.250 - 0.750 GeV have been 
employed for different meson states [7]. Even the case of (3 = m q is adopted 
by Tao Huang et al. in [7]. On the other hand in [10], we have argued that 
the size of the wave packet representing the bound state is important here. 
Therefore we have raised the value of f3 up to 0.6 GeV which we believe is 
quite safe. The behavior of the fragmentation function (7) is shown in Fig. 1 
for (3 = 0, 0.2, 0.4 and 0.6 GeV. We take the maximum value of (3=0.6 GeV 
and use it for our further considerations. 



3 Inclusive production cross section 

We have employed the idea of factorization to evaluate the J/ip production 
cross section at hadron colliders. For pp collisions we may write 
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Fig. 1. Fragmentation function for J/ip production. While the dashed curve repre- 
sents the function when the Fermi motion is off, the solid ones show the behavior of 
this function when the Fermi motion is on. The (3 values are indicated. The function 
picks up as (3 increases. This gives raise to increase in the fragmentation probability. 

1(7 {pp -> c -> J/ip(p T )X) = J2 J dx 1 dx 2 dzf i/p (x 1 , n)fi( p {x 2 , n) 
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Where f\j are parton distribution functions with momentum fractions of X\ 
and X2 (different from x\ and x 2 which appear in (1)), a is the charm quark 
production cross section and D(z,fi,f3) represents the fragmentation of the 
produced heavy quark into cc state with confinement parameter (5 at the 
scale fi. We have used the parameterization due to Martin-Roberts-Stiriling 
(MRS) [11] for parton distribution functions and have included the heavy 
quark production cross section up to the order of a z s [12]. The dependence 
on ii is estimated by choosing the transverse mass of the heavy quark as our 
central choice of scale denned by 



= yPr 2 (parton) 



(10) 



and vary it appropriate to the fragmentation scale of the particle state to be 
considered. This choice of scale, which is of the order of pt (parton), avoids the 
large logarithms in the process of the form ln(mg//x) or ln(pr//i). However, 
we have to sum up the logarithms of order of fiR/rriQ in the fragmentation 
functions. But this can be implemented by evolving them by the Altarelli- 
Parisi equation [13]. The following form of this equation is used here 
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Fig. 2. The differential cross section for direct fragmentation production of J/ip 
and its subsequent decay J ftp — ► ^ ' [i~ at the Tevatron Run I energies. While the 
dashed curve is obtained using (5) or equally (7) with (3 = 0, the solid one is due 
to (7) with j3 = 0.6 GeV. The result is compared with the CDF Run I data. Other 
contributions are not included. The scale is chosen to be 



Here Pq^q(x = z/y,fx) is the Altarelli-Parisi splitting function. The boundary 
condition on the evolution equation (11) is the initial fragmentation function 



[z,n,(3) at some scale \x = ji - In principle this function may be calcu- 



lated perturbatively as a series in a s at this scale. 



Detection of final state requires kinematical cuts of the transverse momen- 
tum, pt, and the rapidity, y. We have imposed the required pip* and y cnt in 
our simulations for different colliders as required and have used the following 
definition of rapidity 
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Fig. 3. The differential cross section for direct fragmentation production of J/tp at 
Tevatron Run II. The two curves are obtained using (7) with (3 = and 0.6 GeV 
respectively. The scale has been set to 2(ir. 

4 Results and discussion 



We have used a light-cone wave function to introduce the Fermi motion in J ftp 
production in direct fragmentation channel and obtained its fragmentation 
function in leading order perturbative regime. In this function the confine- 
ment parameter switches the effect of Fermi motion and within its physically 
acceptable values, the function demonstrates the effect of Fermi motion in an 
interesting manner. The motivation of introducing this effect is to see its role 
in improving the QCD versus experimental results for the J/ip cross section 
at the Tevatron energies. 

In the case of the J/ip state such a fragmentation function, with a reasonable 
choice of the confinement parameter, gives raise to a significant increase in 
the fragmentation probability due to the Fermi motion. Such a behavior is 
illustrated in Fig. 1. 

To see the effect in the production rates, we have used the usual procedure 
of factorization for J/tp production in hadron colliders. We present the pr 
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Fig. 4. The differential cross section for direct fragmentation production of J ftp at 
the RHIC and the CERN LHC. The shift due to the Fermi motion is significantly- 
increased in the case of LHC but shows to be less important at the RHIC. In 
both cases the fragmentation function (7) is employed with (3 = and 0.6 GeV 
respectively. 

distribution of BR(J/ip — > fi + fi~)da/dp T for the cases of the Fermi motion 
off and on along with the CDF data at Run I in the Fig. 2. The branching 
ratio BR(J/ip — > / u + / u _ )=0.0597 is taken from [14]. The poor agreement with 
data is due to the fact that here we have only considered the contribution 
of pp — > c — > J /if) — ■> fi + fi~. Similar behavior at Tevatron Run II energies is 
shown in Fig. 3. We have also extended our study to the cases of the RHIC 
and the CERN LHC pp colliders. Here we provide the p T distributions of the 
differential cross sections for c — > J/ip and compare the two cases of the Fermi 
motion off and on in the Fig. 4. In all cases we have used f3 = 0.6 GeV for the 
confinement parameter in the fragmentation functions. Naturally, the results 
for p in the range of - 0.6 GeV fall between the above results. 

We have also calculated the total integrated cross sections for each case. We 
found that the total cross sections for with and without Fermi motion essen- 
tially remain unchanged within the uncertainties of Mote Carlo simulations. 
The reason is first due to the momentum sum rule which the fragmentation 
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functions should satisfy. In other words although the modification of frag- 
mentation functions by the Fermi motion redistributes the final states, the 
integrated cross sections are left unchanged. Alternatively although the Fermi 
motion increases the fragmentation probability for the state, i.e., introduces 
a state with overall higher mass, the cross section is lowered by just the same 
amount when we introduce the effect in calculation of the total integrated 
cross section. 

First of all we note that the effect of Fermi motion is indeed significant and 
that the implementation of such a study in all production channels of J/if> 
may give raise to considerable enhancement of the color singlet differential 
cross section. 

It is evident from the Figures 2,3 and 4 that the effect increases with increasing 
y/s. The kinematical cuts play important role apart from s/s. The large cross 
section at the RHIC compared with the LHC in the Fig. 4 is an example. 

Finally we consider the uncertainties involved in our study. There are two main 
sources of uncertainties. The first is about the simulation of J/ip production 
at hadron colliders such as the uncertainties along with the fragmentation 
functions and parton distribution functions. These kind of uncertainties are 
well discussed in the literature. The second source of uncertainty is due to the 
choice of the confinement parameter in the fragmentation function. Relying 
on our discussion in section 2, our choice of /3 = 0.6 GeV seems to be justified. 
Future determination of this parameter will shed more light on this situation. 
It is worth mentioning that our choice of charm quark mass, i.e., 1.25 GeV, 
have put our results in their upper side and that the change of the charm 
quark mass in its acceptable range does not have significant impact on the 
Fermi motion effect in J/ifj production. 
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